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ABSTRACT 

Gamma-ray burst time histories often consist of many peaks. These 
peaks tend to be narrower at higher energy. If gamma-ray bursts are cos- 
mological, the energy dependence of gamma-ray burst time scales must be 
understood in order to correct the time scale dependence due to the expan- 
sion of the universe. By using the average autocorrelation function and the 
average pulse width, we show that the narrowing with energy follows, quite 
well, a power law. The power law index is ~ -0.4. This is the first quan- 
titative relationship between temporal and spectral structure in gamma-ray 
bursts. It is unclear what physics causes this relationship. The average au- 
tocorrelation has a universal shape such that one energy range scales linearly 
with time into all other energy ranges. This shape is approximately the sum 
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of two exponential. 
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1. INTRODUCTION 



The Burst and Transient Experiment (BATSE) on the Compton Gamma 
Ray Observatory has deepened the mystery of gamma-ray bursts (GRBs) 
rather than solving it. GRBs appear to be isotropic on the sky yet there 
is a dearth of faint events compared to the brightest events (Meegan et 
al. 1992). The two most likely explanations for this situation are the bursts 
are at cosmological distances (and the dearth of events is due to effects 
associated with the expansion of the universe) or the events are from an 
extended halo about our Galaxy (and the dearth of events is due to a decrease 
in the density of neutron stars in the halo). If cosmological, the expansion 
of the universe shifts the photon energies by a factor of 1/(1 + z) where 
z is the redshift and stretches the temporal structure by a factor of 1 + 
z. Indeed, time dilation has been claimed on all time scales within GRBs 
(Norris et al. 1994, 1995a, Davis et al. 1994). A variety of tests have been 
used to detect the time dilation. The total-count test, wavelet-power test, 
and aligned-peak test consistently require a factor of 2 dilation between the 
dimmest and brightest BATSE events (Norris 1994, Norris et al 1994) as do 
the T 50 and T 90 distributions (Norris et al. 1995a). Norris et al. (1994, 1995a) 
interpret the factor of two dilation as consistent with the GRB Log N-Log 
P distribution, although perhaps with some evolution. Fenimore & Bloom 
(1995) contends that it is not consistent when one includes all the factors 
relating distance to time dilation. One key factor involves the tendency for 
peaks in GRB time histories to be narrower at higher energy. Fishman et 
al. (1992) noted that individual peaks frequently are narrower and better 
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defined at higher energies. Link, Epstein, and Priedhorsky (1993), showed 
that this is a prevalent property of most bursts. In this paper we show 
that there is a well defined relationship for the average width of peaks as a 
function of energy. We, like Link, Epstein, and Priedhorsky (1993) will use 
the autocorrelation function of the GRB time history. We will show that 
the average autocorrelation function for many bursts is a very well behaved 
function with a shape that is universal. Heuristically, an autocorrelation 
measures the average relative intensity between points in the time history 
that are separated by an amount of time called the lag. If a GRB time 
history is stretched by a factor, the width of an autocorrelation function 
that covers the entire time history is increased by the same factor, assuming 
that noise does not dominate. In some cases it is not practical to use the 
entire time history. We always use enough of the time history such that the 
actual length that we use does not have a large impact on the analysis (see 
section 3). As such, it can be used to detect changes in time scales that might 
be associated with the expansion of the universe or measure the average peak 
width as a function of energy. The average autocorrelation is fairly immune 
to systematic effects such as the identification of the highest peak. The 
noise is explicitly accounted for by calculating the expected autocorrelation 
given the noise level. The average autocorrelation is similar to the aligned 
peak tests in that the peak of each burst is used as a fiducial to form an 
average and is sensitive to time scales the order of a few seconds. However, 
the autocorrelation uses much more data so it has better statistics. It is 
roughly equivalent to aligning most of the peaks in a burst rather than just 



4 



the highest. 

2. INSTRUMENTATION 

The BATSE experiment on CGRO uses eight large area detectors (LADs) 
to locate and study GRBs over a large dynamic range. Usually, the detector 
count rate is recorded with modest temporal resolution (1.024 sec). When 
there is a statistically significant increase in the counting rate above the back- 
ground, special modes are triggered and the data is recorded at a variety of 
temporal and energy resolutions. See Fishman et al. (1992) for further de- 
tails on the instrumentation and data modes. For the purposes of this study, 
we will use the four channel triggered data. This consists of the counts in 
four broad energy bins labelled "1" for 25 to 57 keV, "2" for 57 to 115 keV, 
"3" for 115 to 320 keV, and "4" for above 320 keV which is effectively 320 
to 1000 keV. A data set labelled "1+2" combines "1" and "2" together to 
effectively create a 25 to 115 keV channel. A memory records ~ 2 sec of data 
before the trigger and the duration of the recorded data after the trigger is 
~ 240 sec. The time resolution for this data is 0.064 sec. In Norris et al. 
(1994), the period prior to the trigger is augmented by rebinning the contin- 
uously available 1.024 sec samples into 0.064 sec samples. This extends the 
pretrigger by ~ 16 sec. 

We will use the same data set as used by Norris et al. 1994 including 
the augmented pre-trigger. The on-board electronics adds data from all the 
detectors that triggered, this can vary from two to four detectors. Bursts 
were assigned by Norris et al. (1994) to a brightness class based on the 
largest net count rate in 0.064 sec samples in channel 1+2+3+4. Those 
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bursts with count rates between 18,000 and 250,000 cts sec - are called 
"bright" bursts, those with counts rates between 2,400 and 4,500 cts sec -1 
are "dim", and those with count rates between 1,400 and 2,400 are called 
"dimmest". (Events with an intermediate count rate are not used since the 
time dilation effects are largest for well separated classes.) Short events 
(defined here and in Norris et al. 1994 to have durations less than 1.5 sec) 
were excluded from the study. In this paper, we seek the intrinsic variation 
with energy of the width of the temporal peaks. We use only the bright events 
to avoid potential effects due to different distances including time dilation 
from the expansion of the universe. Even if GRBs come from cosmological 
distances, under the standard candle assumption these events are all from 
approximately the same distance and therefore have the same stretching due 
to the expansion of the universe. There were 45 useable events in the bright 
class. This data set is the same as used by Norris et al. 1994 although 
processed completely independently with the exception of the selection of 
events and the augmentation of the pre-trigger data. 

3. THE AUTO-CORRELATION FUNCTION 

The autocorrelation function for GRBs was investigated by Link, Ep- 
stein, and Priedhorsky (1993) where it was shown that time scales are almost 
always shorter at higher energies. Following Link, Epstein, and Priedhorsky, 
let m, be the observed gross counts in discretely sampled data in n bins of 
equal size AT ranging from — nAT/2 to +nAT/2, about the largest peak in 
the GRB time history. Here rrii is number of counts, so it follows Poisson 
statistics. Let hi be the corresponding background counts. We determined 
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the background by a linear fit to regions before and after the bursts where 
it was judged by eye to be inactive. The net counts are Cj = m, — bi and the 
estimate of the true autocorrelation as a function of the lag, r = j AT, is 



A(r) = 1 



3 = 



n/2 



n/2 < j < n/2 



(1) 



i=-n/2 



A 



Here, c^ + j is zero if i + j > n/2 or i + j < —n/2. When studying individual 



events (as Link, Epstein, & Priedhorsky [1993] did), one can always use the 
entire time history. We average many events together and must use them 
in a uniform manner. This requires us to select the duration to use. If we 
selected a very short duration then the auto-correlation is not sensitive to 
time stretching (a flat-topped burst exceeding the selected duration would 
show no difference in the auto-correlation). In we selected a very large du- 
ration, then there would be many instances when bursts would have to be 
left out of the averaging because the peak occurs within ±AT/2 such that 
Ci+j is not defined for part of the needed range. Although bursts can be very 
long, usually emission more than a few seconds away from the largest peak 
contributes only a little to the auto-correlation function. We have found only 
small differences for all nAT > 8 s and have used nAt = 16 s throughout 
this paper. By definition, the autocorrelation is symmetric, A(r) = A(—t). 
The normalization, Aq, is 



n/2 



- m % 




i=-n/2 
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An autocorrelation without normalization would have a large peak at 
r = where all the noise adds coherently and would be count rate dependent 
at r 7^ 0. The -m, term in equation (2) normalizes the autocorrelation to 
that expected without noise. 

In order to fit functions to the observed autocorrelations, we require 
a measure of its uncertainty. The terms of A(t) are not statistically inde- 
pendent. However, our use of the uncertainty is only to obtain a relative 
goodness of fit. In fact, each term of the variance on A{r) will be approx- 
imately the same so its exact value is not important. The variance of the 
numerator of equation (1) is 

n/2 

Oc-kCj = C i \ C i+j\ + l C *l C i+j • (3) 

i=—n/2 

(We assume that the variance propagated from the background is small since 
the background is based on much more data that the individual points.) The 
variance on the normalization is 

n/2 

°\= E 4c '+™, • ( 4 ) 

i=-n/2 

Combining equations (3) and (4) gives the variance on the j-th term of the 
autocorrelation: 

a A{ 3 AT) ~ A 2 + A 2 °A ■ \p) 

The average of a fair number of GRB autocorrelations is quite stable 
and shows only a small variation. Let A(i,r) be the average autocorrela- 
tion for the 2-th channel or combination of channels. Figure 1 shows the 
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average of the bright events for the four channels of the LAD data, that is, 
A(i,r) = Ylk=i Ak(h T )/NB where Nb is the number of bright events (45) 
and k denotes different bursts. The normalization of each autocorrelation 
is such that each burst contributes equally to the average autocorrelation 
independent of its brightness. Note how clear the energy dependency is in 
Figure 1. Figure 1 is semi-log and the curves appear nearly as straight lines 
so the shape of the autocorrelation is approximately an exponential. 

Each energy channel is nearly an exact time stretched version of the 
others. We define Sij to be the best-fit factor that scales A (i, r) into A (j, r). 
It is found by minimizing 



Here mAT is the range of lags that is used in the fit. This range is set by 
where the functions are well defined. The auto-correlation of the highest 
energy channel begins to have significant noise at ~ 2.5 s (see Fig. 1) so we 
have used mAT = 2.5 sec. Since A(t) is symmetric nothing is gained by 
including negative lags. 

In each panel of Figure 2, we show the average autocorrelation for two 
channels. For example, in Figure 2a, we show A(1,t) (i.e., 25 to 57 keV) 
and A(2,t) (i.e., 57 to 115 keV). Also shown as a bold curve is the time 
stretched autocorrelation of the higher energy channel that best fits the lower 
energy channel (e.g., the bold curve in Fig. 2a is X21A (2, £2,1 t)). The over- 
all scaling (A) is always very near unity. Note, for example, in Figure 2c 
that the bold curve slightly exceeds unity at r = 0. For energy channels 
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with poorer statistics, the uncertainty of the normalization is reflected in cr 2 _ 

A 

which is nearly constant as a function of r since mAT (2.5 sec) is much less 
than nAT (16 sec). The goodness-of-fit parameter in equation (6) will not 
follow the x 2 statistic because the points are not independent. The purpose 

of cr 2 _ is to balance the uncertainty in the overall scale factor (A) with the 

A 

uncertainty in the time stretching (S). Figure 2 demonstrates that the av- 
erage autocorrelation has a universal shape (but different time stretching) 
for all energies. Note how well the time-stretched higher energy autocorre- 
lation always agrees with the broader (lower energy) autocorrelation. Even 
the highest energy range (320-1000 keV), which showed a deviation from an 
exponential in Figure 1 scales exactly into the lower energy autocorrelations. 

The curves in Figure 1 are not pure exponentials, there is a slight curve 
to the histograms. We have tried to fit a variety of functions to the his- 
tograms. A single exponential (exp -T / T °) fits poorly, especially the higher 
energy channels. A function such as exp aT + fits the lower energy channels 
well but not the higher energy channels. Although not unique, the most 
successful function that we tried is a sum of two exponentials: 

A(i, t) = fii exp-^l/"* 1 +(1 - fii) exp-l T l/ Q ' 2 (7) 

where i denotes the 4 energy channels. To determine the free parameters in 
equation (7), we minimize 



n/2 r n2 

* 2 = E 



A(i,jAT) -ftexp-l 1 "!/*' 1 -(1 - Pi) exp-l T l/ Q * s 



(8) 



The parameters an and ot.ii are found by searching the parameter space and 
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Pi is found analytically from <5% 2 /Sfii = 0. The best fits are 

A(1,t) = 0.66exp-l T l/ 2 - 40 +0.34exp-l T l/ 25 , (9a) 

A(2,t) = 0.64exp-l T l/ 1 - 74 +0.36exp-l T l/ 25 , (96) 

A(S, r) = 0.48 exp-l T l/°- 94 +0.52 exp"^/ 9 ' 9 , (9c) 

A(4,r) = 0.53exp-l T l/°' 56 +0.47exp-l T l/ 6 ' 5 . (9d) 



All values of above 25 are equally consistent with the data. These curves 
are plotted as curves in Figure 1. Equation (8) does not follow the % 2 statistic 
because the A terms are not independent so we cannot qualitatively evaluate 
the fit. However, as seen in Figure 1, the fit is excellent. 

4. ENERGY DEPENDENCE OF TIME SCALE 
We will characterize the energy dependence of the typical time scale in 
the GRB time history using two different measures: the width of A and 
the width of the average pulse profile (from Norris et al. 1994, 1995b). The 
solid triangles in Figure 3 are the width (W ac ) of each autocorrelation from 
Figure 1 as measured by where \thA(t) = 0.5. Since Figure 3 is log-log and 
the points nearly lie on a straight line, we have fit a power law to the points. 
The best fit power law is 

W ac (E) = 17.4£-°' 43 . (10) 

This function is shown in Figure 3 as a solid line. This is a robust result. 
Using the width at other values of A gives similar results. Also, the fact that 
the auto-correlation function for each energy can be scaled into another and 

11 



they overlap so well (Fig. 2) implies that the power law holds for more than 
just the point where h\A(r) = 0.5 (see discussion of eq. [14]). 

One thing that is not clear in our formulation is what energy to place 
the points at. We have placed them at the energy corresponding to the lower 
energy bound of the channel they represent. The autocorrelation function 
is quadratic in counts (see eq. [1]) so for any particular channel, the width 
reflects where most of the counts are. For example, if one combines channel 
3 and 4, it has effectively the same width as 3 only. If we were to use the 
midpoint of the channel, it is still a power law: W ac = 18.1E~ 0A0 . Another 
possible measure of the time scale of A is how much one energy range needs 
to be stretched to map it into another energy range. This is not independent 
from W ac but serves as another way to measure it (eq. [6] rather than eq. [1]). 
Using channel 1 as a baseline (i.e., if S\\ = 1), S 1 ^ 1 = 0.78, = 0.54, and 
S^ 1 = 0.33 give how much the autocorrelations of the higher energy channels 
are narrower as a function of energy Fitting the ^.j^ 1 points gives 



A second measure of the time scale of GRBs comes from the average pulse 
width. Norris et al. 1994 decomposed GRB time histories into individual 
pulses and found the average rise and fall time scales. The width (in seconds) 
of the rise/fall of the average pulse profile are 0.22/0.44, 0.17/0.32, 0.13/0.27, 
and 0.08/0.18 for BATSE channels 1, 2, 3, and 4, respectively. We plot in 
Figure 3 as squares the sum of the rise and fall times as a function of energy. 
Again, the energy dependence of the time scale appears to be a power law. 




4.45£: 



-0.46 
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In this case, the average pulse width is 




-0.37 



(12) 



which is plotted in Figure 3 as a dashed line. Another measure of the pulse 
width is the average full width, half maximum (FWHM). Norris et al. (1995b) 
reports that the average FWHM for the four BATSE channels are 0.817, 
0.616, 0.473, and 0.287 sec. These widths can be fit by a power law as well: 



This, too, is very robust. Norris et al. (1995b) reports the width for seven 
different fractions of the peak height and all seven can be fit with a power 
law. 

We note that the average pulse width is much less than the width from 
the autocorrelation function. The individual pulses are narrower than the 
clusters of peaks that often determine the autocorrelation width. However, 
there is not a simple relationship between the two measures. For example, 
simulations of shot noise with pulses the order of the average pulse width 
produce autocorrelation functions that are much narrower than observed. 

In equation (8), we fit each energy channel separately using 12 param- 
eters. From Figure 3, we see that the energy dependence of the time scale 
in GRBs is a power law. The parameters an and a^ found in equation (9) 
do not follow a power law. However, it is possible to have a functional form 
that has a power law dependency on energy and fits within the noise. We fit 
all four curves in Figure 1 with 



Wfwhm = 



-0.42 



(13) 



A (i, r) = (3 exp 




+ (1 - (3) exp 
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This form accommodates the three characteristics of the average autocor- 
relation: it consists of two exponentials, the energy dependence scales as 
a power law, and the shape of one energy range scales linearly with time 
into all the other energy ranges. In Figure 4 we show the autocorrelations 
from Figure 1 fit with equation (14). The best fit parameters are (3 = 0.55, 
k\ = 8.75, &2 = 154, and a = 0.45. Although the curves deviate some from 
the histograms, the fits are within the uncertainty expected for the average 
autocorrelation caused by variations in the choice of bursts that are included 
in the average. Thus, we consider the fit acceptable. 

In summary, we find that the average autocorrelation of GRBs time 
histories is a universal function that can measure the time scale as a func- 
tion of energy. The dependence is a power law in energy with an index 
that is between 0.37 and 0.46, depending on how it is measured. This is 
the first quantitative relationship between temporal and spectral structure 
in gamma-ray bursts. The energy dependence is important for two reasons. 
First, the shape may indicate the underlying physics responsible for the time 
history. For example, the subpeak's temporal width might be produced by 
the growth of a shock within a relativistic expanding shell in a cosmolog- 
ical GRB and the power law dependence on energy is related to how the 
shock converts bulk motion into gamma-rays. Alternatively, the energy de- 
pendence might be related to how a disturbance propagates on the surface of 
a neutron star. Second, in order to interpret the time dilation due to the ex- 
pansion of the universe, one must understand the energy dependence which 
competes with the cosmological time dilation to form the temporal width. 
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Fenimore & Bloom (1995) includes the energy dependence in the interpreta- 
tion of time stretching as a function of burst brightness and concludes that 
the observed time dilation is not consistent with the observed log TV-log P 
distribution (Fenimore et al. 1993) unless there is strong evolution and it is 
only coincidental that the log TV-log P distribution shows a -3/2 power law.. 
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FIGURE CAPTIONS 



Fig. 1 The average autocorrelation of 45 bright BATSE gamma-ray bursts 
in four energy channels. At higher energies, gamma-ray bursts have shorter 
time scales. The solid curves are fits of the sum of two exponentials to the 
autocorrelation histograms. 

Fig. 2: Comparisons of pairs of energy channels for the BATSE bright 
events. Each panel shows the average autocorrelation for two energy channels 
from the BATSE LAD data. The bold histogram is a best fit, time stretched 
version of the narrower (higher energy) channel fit to the broader (lower 
energy) channel. The average autocorrelation apparently has a universal 
shape which is approximately exponential. 

Fig. 3: The energy dependency of the time scale as determined by the 
average autocorrelation and the average pulse width. The triangles are the 
half width of the autocorrelation function and the solid curve is a best fit 
to the triangles. The squares are the sum of the rise and fall of the average 
pulse profile and the dashed line is a best fit to it. In each case, the time 
scale of the temporal structure within the GRB scales as a power law of the 
energy with an index of ~ -0.4. 
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Fig. 4 The average autocorrelation of 45 bright BATSE gamma-ray bursts 
in four energy channels fit with the sum of two exponentials where widths 
scale as a power law. 
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Fig. 1: Energy Dependency of 13ATSE Bright 
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